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ABSTRACT 

At high energy densities any quantum field theory is expected to have an effective hydrodynamic 
description. When combined with the gravity/gauge duahty an unified picture emerges, where 
gravity itself can have a formal holographic hydrodynamic description. This provides a powerful tool 
to study black holes in a hydrodynamic setup. We study the stability of plasma balls, holographic 
duals of Scherck-Schwarz (SS) AdS black holes. We find that rotating plasma balls are unstable 
against m-lobed perturbations for rotation rates higher than a critical value. This unstable mode 
signals a bifurcation to a new branch of non-axisymmetric stationary solutions which resemble a 
"peanut-like" rotating plasma. The gravitational dual of the rotating plasma ball must then be 
unstable and possibly decay to a non-axisymmetric long-lived SS AdS black hole. This instability 
provides therefore a mechanism that bounds the rotation of SS black holes. Our results are strictly 
valid for the SS AdS gravity theory dual to a SS gauge theory. The latter is particularly important 
because it shares common features with QCD, namely it is non-conformal, non-supersymmetric 
and has a confinement/deconfinement phase transition. We focus our analysis in 3-dimensional 
plasmas dual to SS AdSs black holes, but many of our results should extend to higher dimensions 
and to other gauge theory /gravity dualities with confined/deconfined phases and admitting a fluid 
description. 
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1 Introduction 



The equations of fluid dynamics and properties of fluids at large have been used for centuries, 
not only to describe fluids but also as analogue models for other more complex phenomena. For 
instance, early experiments with liquid drops by Plateau [1] were aimed at understanding the effect 
of gravity on planets (surface tension was then a model for the gravitational force). Another well- 
known example, is Bohr and Wheeler's p| proposal to describe nuclear fission as the rupture of 
a charged liquid drop, where now the surface tension plays the role of nuclear forces. In general 
relativity, the membrane paradigm, whereby a black hole horizon is mimicked by a stretched fluid 
membrane, provides another example of the power of analogue models, with useful applications in 
astrophysical systems. Still in a gravity setup, it was recently suggested to use fluid analogs to 
explain phenomena observed in general relativistic scenarios, in particular the classical instability 
of black strings and branes [3] . Accordingly, the gravitational Gregory-Laflamme instability would 
have a counterpart in the Rayleigh-Plateau instability of fluid mechanics [H [5] (responsible for the 
breakup of liquid jets and tubes). 

1.1 Dual hydrodynamic description of gravity 

The anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence adds an interesting twist 
to this story, making these analogies powerful and formal. Indeed, it has emerged with the work of 
Fermi [6], Landau [7j and others, that often the complicated time dependent dynamics of quantum 
fields is approximated by a fluid model description. It seems in fact that any gauge theory has 
a hydrodynamic limit [8, 9j. Combining these ideas together one expects gravity to have a dual 
hydrodynamic description [9j. This expectation has been formally verified in \10\ lllj . and later 
in [121 I13j . where it was explicitly shown that a gravitational geometry satisfies perturbatively 
Einstein-AdS gravity to any order as long as the associated holographic stress tensor T^'^ (read 
from the AdS/CFT dictionary) has vanishing spacetime divergence, ^^T^^^ = 0. One recognizes 
the latter equations as those that govern fiuid dynamics. At leading order the stress tensor is 
that of a perfect fluid; in the next-to-leading order in the perturbation, T^^ gets a contribution 
that describes viscosity and dissipation effects; at higher order T^^^ provides information about the 
fluid relaxation timescales. Thus Einstein-AdS gravity is indeed dual to hydrodynamics, in the 
appropriate regime. 

The hydrodynamic description of gravity (THl [HI [ISl Ej has support on the AdS / CFT duality, 
relating type IIB string theory on AdSs x with = 4 Super Yang-Mills (SYM) gauge theory. 
SYM differs considerably from QCD. For example, as opposed to SYM, QCD is non-conformal, 
non-super symmetric (non-SUSY) and has both a conflned and a deconfined phase. Thus, the 
holy graal of the field is to find a string/QCD duality, which would allow one to study hard non- 
perturbative phenomena in QCD through a weak-coupling perturbative analysis of the dual string 
system, and vice- versa. So far this programme has not been completed (see |141ll5j for discussions), 
however some gravity/gauge dualities are known where the gauge theory shares some important 
common features with QCD. The simplest example is the Scherk-Schwarz (SS) compactification 
of a 5-dimensional (5(i) CFT which yields a Ad non-conformal, non-SUSY gauge theory with a 
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confinement /deconfinement phase transition [Jiip. The original CFT is 5d maximally SUSY SYM 
theory that describes field excitations living in a stack of Z)4-branes. Identifying periodically one of 
the worldvolume directions of the L>4-branes, imposing anti-periodic boundary conditions for the 
fermions along this direction, and finally dimensionally reducing along this compact direction one 
gets the desired SS gauge theory. The gravitational dual description of this system is obtained by 
taking the appropriate decoupling limit of the geometry describing the near-extremal D4-branes 
with the compact SS worldvolume direction. Because of the presence of this SS direction there 
are two solutions: one is a black brane (deconfined gluon phase) and the other one is the AdS- 
soliton (confined hadronic phase) [l9j. The solution that dominates the partition function is the 
one that minimizes the free energy of the system. One finds that there is a critical temperature 
Tc - the confinement temperature - above (below) which the black brane (AdS-soliton) minimizes 
the free energy [16]. The confinement temperature is the one where the Euclidean time circle 
has the same length as the SS circle. At this confinement temperature the two phases can co- 
habit in equilibrium separated by a domain wall. Not less important, at this temperature one 
can have a confinement /deconfinement phase transition in the gauge theory which corresponds in 
the dual gravity side to a phase transition between the thermal AdS-soliton and the black brane 
phases |16j (in the context of global AdS backgrounds the transition between thermal global AdS 
and the Schwarzchild black hole solution is known as Hawking-Page phase transition). Some of 
these properties are schematically represented in Fig. [TJa). The black hole, solution of the SS 
compactification of AdSs gravity on a circle, that interpolates between the black brane phase and 
the AdS-soliton confined phase with a domain wall in between is still not known. But in Ref. [20] 
a numerical solution was found which describes an infinite planar domain wall separating the black 
brane in one side from the confined AdS-soliton on the other, at the confinement temperature. This 
solution is expected to describe approximately the near-horizon geometry of the above mentioned 
black hole solution in the limit where the black hole is large. 

Strong arguments suggest that at high energy densities the SS compactification of 5d CFT 
also has a long wavelength effective hydrodynamical description [20j. Indeed, as described two 
paragraphs above this is certainly true for a CFT, and a similar proof (although necessarily more 
complicated) should follow similarly for a SS compactification of a CFT. Then, the corresponding 
gravity/gauge duality, asserts that 3d fiuid dynamics provides an effective theory describing the SS 
compactification of AdS^ gravity in the long wavelength regime [201 121j . That is, on the boundary 
of a SS compactification of AdSs (asymptoting to Ai^ x S^, with the distinguished SS circle), 
the black hole is described by a plasma lump immersed on the vacuum confined phase with a 
domain wall with surface tension separating the two phases. This is schematically represented in 
Fig. dlb). In this description, the black hole horizon maps to the full plasma lump bulk, not to its 
boundary. The SS circle plays a minor role on the fluid description (meaning that the plasma lumps 
are translationally invariant along this direction) but as we enter through the radial holographic 
direction into the bulk, the SS circle must shrink to zero size at the horizon where the domain 

^SS theory has a transition from a hadronic phase to a gluon phase. To have instead a quark-gluon phase at 
high temperature we need to add fundamental matter to the model. This can be done through the introduction of 
probe D6-branes [17] or D8-branes [18]. The latter system is known as holographic QCD or Sakai-Sugimoto model. 
Currently, these gauge systems are the closest to QCD we can have with a theory that has a gravity dual [15] . 
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wall is. This implies that the topology of the corresponding event horizon is given by the fibration 
of the SS circle over the plasma lump geometry, with the circle shrinking to a point on the 
boundary. So, e.g., a plasma ball with a disk topology corresponds in the bulk to a black hole 
with horizon topology S^, and a plasma ring with topology S"^ x / is dual to a black ring with 
topology X S"^ [21J. To leading order, i.e., without dissipation, the plasma lump is described 
by a perfect fluid stress tensor with an equation of state characterizing the fluid from which the 
gauge theory is "made of". The domain wall contributes with a boundary term, proportional to 
its surface tension, to the stress tensor. Finally, the system is calibrated in such a way that the 
confined phase exterior to the plasma lump is vacuum with zero pressure. These plasma balls and 
plasma rings were studied in great detail in [21] with an emphasis on the AdSs case, and the full 
phase diagram for balls and rings in AdSg was obtained more recently in (22j . 
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Figure 1: (a) Dominating phases in the gauge theory and in gravity. For temperatures below the 
confinement temperature Tc, the confined hadronic phase (AdS-soliton in gravity) dominates the 
partition function, while above Tc the deconfined gluon phase (black brane in gravity) dominates. 
At Tc the system suffers a first order confinement/deconfinement phase transition, (b) At Tc and 
in its vicinity, in the dual fluid description on the boundary of the AdS gravity solution, the black 
hole is described by a plasma phase immersed in the confined vacuum phase and separated by a 
domain wall with surface tension. 



1.2 A brief summary of our results 

In the present study we analyze the stability of these plasma lumps, with an emphasis on plasma 
balls. We find that rotating plasma balls are unstable against m-lobed perturbations (with m 
being the azimuthal number of the perturbation) if their rotation is higher than a critical value. 
We further find that the marginal unstable mode is a bifurcation point to a new branch of stationary 
solutions in the phase diagram of solutions. This new phase describes m-lobed plasma lumps. In 
the simplest m = 2 case we have a 2-lobed configuration which presumably (if we draw from 
experience with classical incompressible fluids [23 ^ 124 ^ [25]) goes over to a peanut-like configuration 
for large enough angular momentum. In this work, we shall refer to such solutions interchangeably 
as rotating plasma peanut, or 2-lobed configurations. A phase diagram including also this new 
family is sketched in Fig. [2l The associated gauge/gravity duality will then be used to predict 
that black holes asymptoting to a SS compactification of AdSs should also be unstable to m-lobed 
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perturbations. The result of the current study provides a good example of how the hydrodynamic 
description of gravity can provide a powerful predictive tool to discuss black hole physics and the 
associated dual gauge theory. Note that although we restrict our analysis to d = 3 plasmas dual 
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Figure 2: Schematic phase diagram with angular velocity vs angular momentum of stationary 
solutions at fixed energy. The rotating plasma balls and plasma rings ( "fat" and "thin" rings merge 
at point M) phases were already studied in Ref. [21]. Here we find that above a critical rotation 
(represented by points A, B, ...), plasma balls become unstable against m-lobed perturbations (the 
unstable balls are represented by a dash line). The marginal point where the instability sets in is 
a bifurcation point to a new branch of non-axisymmetric m-lobed plasma lumps. In the simplest 
m = 2 case, this is a rotating plasma peanut. This phase diagram summarizes some of our main 
results, accurately represented in Fig. [71 



to black objects in SS AdSs, our main results and conclusions should extend to any d-dimensional 
plasma lump dual to SS AdS(i+2 black objects, and to other gauge/gravity dualities with confined 
and deconfined phases and admitting a fluid description (see however discussion in Section [6]). 

The instability of rotating plasma balls could have been guessed from classical works and ex- 
periments with rotating fluids [U [231 [231 [26[ [27] . -A- similar reasoning was recently used [131 [28] to 
relate the Gregory-Laflamme gravitational instability of black strings and branes to the Rayleigh- 
Plateau instability of plasma tubes. In the rest of this introduction we describe with some detail 
the aforementioned classical studies on rotating fluids. 

The stability of incompressible Ad non-relativistic fluids was studied in detail by Chandrasekhar 
|23j and later completed by Brown and Scriven [24j, for fluids held by surface tension and in 
the absence of gravity. The analysis was extended in Ref. [5] to a general number of spacetime 
dimensions. In short these works found the following. Start with a static fluid ball and slowly add 
rotation. The fluid ball starts to flatten at the poles. One finds a critical rotation frequency f^c, after 
which two or more families can co-exist. One of the families is the axisymmetric one. Proceeding 
along this family while increasing rotation, one finds a second critical frequency $7* at which the ball 
"pinches" at the origin, and becomes a doughnut like configuration. The axisymmetric family was 
found to be unstable to small perturbations after the point f^c- The other family that bifurcates off 
at Qc is a "two-lobed" configuration (see Figs. 1 and 2 in [5]). Close enough to the bifurcation point, 
this family is stable to small perturbations. There are other families, three-lobed, four-lobed, etc. 
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branching off at different points as one increases rotation. These seem to be always unstable. Recent 
accurate experiments |26j have confirmed the existence, bifurcation points, and stability properties 
of these families. This kind of evolution diagram was found to hold also in gravitationally bound 
(Newtonian) objects with the same qualitative behavior: the Mac-Laurin configuration being the 
axisymmetric family and the Jacobi sequence, a tri-axial ellipsoid, branching off at the bifurcation 
point [291 

The plan of the rest of the paper is the following. Section[2]reviews the relativistic hydrodynamic 
equations, the equilibrium conditions, equation of state, and the conserved charges, that govern 
a Scherk-Schwarz plasma in a 3d Minkowski background (we try to be self-contained) . Section [3] 
then discusses the properties of the axisymmetric lumps of the theory, namely the plasma balls and 
plasma rings, and discusses the regime of validity of the hydrodynamic description. In SectionlHwe 
discuss in detail the stability of plasma balls, the stability of plasma rings is deferred to AppendixlAl 
In Section [5] we find that the marginal stable mode found in the previous section is a bifurcation 
point to a new branch of lobed lumps. Section [6] discusses the consequences of these findings to the 
phase diagram of black hole solutions in SS AdS gravity. 

2 Relativistic hydrodynamic equations 

Here we review the relativistic hydrodynamic equations governing a Scherk-Schwarz plasma in a 3d 
Minkowski background. We also derive the dissipation contribution to the hydrodynamic equations 
which are not easy to find in the literature. We will be interested in plasma configurations in 
mechanical and thermodynamic equilibrium. We follow closely |2H I28j. 

2.1 Relativistic hydrodynamic equations 

Fluid mechanics is an effective description at long distances, valid when the fluid variables vary on 
distance scales that are large compared to the mean free path l^fp of the system. As a consequence 
it is natural to expand the stress tensor in powers of derivatives of the four- velocity u'^. To zeroth 
order in the derivative expansion, Lorentz invariance and the correct static limit uniquely determine 
that the stress tensor is the sum of the perfect fluid plus boundary contributions, 

T^f = {P + P) ^'u"" + Pr , = -<Thndf\ 5{f) . (2.1) 

Here, is the fluid velocity, p, P and a are the density, pressure and surface tension of the fluid. 
The fluid boundary is defined by f{x^) = 0, it has unit spacelike normal = d^f /\df\, and 
= g^'^ — n^^n" is the projector onto the boundary. 

The first subleading order of the derivative expansion introduces dissipation effects in the prob- 
lem. Lorentz invariance and the physical requirement that entropy variation is non-negative de- 
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mands that the dissipation stress tensor is [30j 

T^.l = -Ci?P^^ - 2r^a'''' + qf^u" + u^^q" , (2.2) 
where C is the bulk viscosity, rj is the shear viscosity, k is the thermal conductivity, and (for d = 3) 

a^"' = ^ (p^^Vau'^ + P^^Vau'^) - ^^?P''^ = -KP^'^(a^T + a^T) , (2.3) 

are the acceleration, expansion, shear viscosity tensor, and heat flux, respectively. The last equation 
is the relativistic Fourier law. 

The hydrodynamic equations describe the conservation of the stress tensor, 

^AtZ, + tZ. + tZ^=^- (2.4) 



M \^ pcrf bdry diss j 

It then follows that in the presence of dissipation, the relativistic continuity, Navier-Stokes and 
Young-Laplace equations, are respectively given bjo 

n^V^p +{p + P)d = - q^'a^ - V^g^ + 2ry(j'^'^V^n, , (2.5) 
{p + P)a- = -PA'-V^P + C (P'^'V^^ + iJ^z^V^n'^) + 27? (V^a^^ - nV^"V^tx«) 

- {qf'V^u'' + ^q" + u'^Vf.q" - qf'a^u") , with P'"' = gf"" + uf^u" , (2.6) 



< 

= aK, with K = h^V^n^', (2.7) 



where P^'^ is the projector onto the hypersurface orthogonal to u^, K is the boundary's extrinsic 
curvature, and [Q]> = <5< — <5> is the jump on a quantity Q when we cross the boundary from the 
interior into the exterior of the plasma (we will be interested in the case where the plasma object is 
immersed in vacuum; then the outside contribution in the Ihs of (j2.7p vanishes). In the derivation 
of Eq. (j2.7p . the constraint that the fluid velocity must be orthogonal to the boundary normal is 
used (this guarantees that the fluid is confined inside the boundary), 

u^n^ = 0. (2.8) 

For a conformal plasma it is well-known that the bulk viscosity coefficient vanishes, = 0, and 
that the ratio of the shear viscosity r] to the entropy density s of the plasma is ^ = A-. However, 
our system is described by a Scherk-Schwarz plasma that is not conformal and the dissipation 
coefficients satisfy the relations > 0, ^ > A. and k > 0. 



At leading order the entropy current density is (J5)pcrf ~ su'^ and is conserved. In the first subleading order, 
one gets the extra dissipative contribution (Jg)diss ~ The entropy current density Jg = (J^)perf + (Jg)diss is no 
longer conserved and satisfies TV^ Jg = -^-^^ + C^^ + 2riG^v(y^" > as long as 77, and k axe positive parameters, as 
we assume [2T]. In equilibrium, V^j Jg must vanish. It follows that, q**, 9 and a*"^ must vanish in equilibrium. 

^To go from (|2.4p into (|2.5p - (|2.7[) note that (|2.4p reduces to a volume and a boundary contributions. The latter 
gives the Young-Laplace equation while the former yields the Continuity and Navier-Stokes equations. To get these 
we use V^e(-/) = - (9^,/) ,V ^,5{f) = ~\df\n^,5{f), u^v!' = -1, /i'"'nf, = 0, P^'u^ = 0, u^o!' = 0, u'^n^, = 0, 
q*^Up. — 0, a'"' Up — 0, M''itp = —1, n'^n^ = 1, and V aO'^'^ = 0. We also use properties of the type: n^ny = 1 implies 
n^nf V^n" = 0. 
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2.2 Hydrodynamical and thermal equilibrium conditions. Equation of state 

In this subsection we briefly review some general results derived in Ref. [28] valid for equilibrium 
and general (non-)axisymmetric plasma configurations. We focus our discussion on fluids in a 3d 
Minkowski backgrounqj with stationarity timelike Killing vector ^ = dt and spacelike Killing vector 
X = d<t>i but the results of |28] were derived for a general background geometry. 

A fluid with local entropy density s and local temperature T satisfies the Euler relation and 
the Gibbs-Duhem relation given, respectively, by 

p + P = rs, dP = sdT , (2.9) 

where the latter relation follows from differentiating the former relation and use of the first law 
of thermodynamics. Using these relations and demanding hydrodynamic (i.e., mechanical) equi- 
librium we find that the plasma must also be in thermodynamic equilibrium [28]. So, equilibrium 
plasma configurations must satisfy the hydrodynamic equations discussed in the previous subsec- 
tion with vanishing subleading dissipation and diffusion contributions. It then follows that any 
stationary fluid configuration with local temperature T must have a velocity given by 

u = '^i^-nx). (2.10) 

with constant T and fi. So stationary configurations are rigidly rotating equilibrium solutions with 
constant plasma temperature T related to the local temperature T by the Lorentz factor (we use 
= -1), 

which is the redshift factor relating measurements done in the laboratory and comoving frames. 

Combining the Euler relation (j2.9p and the Young-Laplace equation (equation (j2.7p without the 
dissipative terms), we can relate the plasma temperature T to a combination of several magnitudes 
at the fluid surface, 

r="^ + ^ (2.12) 

We see that T is not simply proportional to the surface tension or to the mean curvature, although 
it grows linearly with both [2^. For a static fluid K will be constant over the surface, but in a 
stationary configuration K is not a constant along the boundary. In the duality to a black hole, T 
corresponds to the Hawking temperature of the horizon. 

A universal behavior of fluids is that they always pick boundary configurations which reduce 
their potential energy for a fixed volume. For static solutions, this implies that the area of the fluid 
surface is minimized. For stationary solutions, the potential energy not only has a surface tension 
term but also a centrifugal contribution. In Ref. ^28j it was shown that this variational principle 
still holds for relativistic fluids and that the minimization problem is equivalent to maximize the 
plasma entropy while keeping its energy and angular momenta fixed. In the gravitational dual 
system black holes satisfy the variational principle that their entropy is extremized for fixed energy 

*We use polar coordinates (t,r, 0), and the non- vanishing afiine connections are F^^ = —r and Vf^ = F^^ — 1/r. 
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and angular momenta. In the duality between SS-AdS black holes and fluid lumps, the entropy, 
energy and angular momentum are identified on both sides, while the temperature is mapped 
according to ()2.12p . Note that unexpectedly (because it looks a priori to be a contradiction), the 
analysis of Ref. [28] shows that maximization of the black hole horizon area is equivalent, for static 
configurations, to minimization of the fluid surface area. The reason for this equivalence can be 
traced back to the fact that the black hole horizon is mapped to the entire volume of the plasma 
and not to the plasma boundary. 

The results quoted so far are independent of the equation of state for the fluid. The gravita- 
tional/hydrodynamic duality where we frame our analysis requires however a particular equation of 
state. Indeed, we are interested in the long wavelength limit of a Scherk-Schwarz compactification 
of a 4d CFT. The 3d (non-conformal) plasma that results from the dimensional reduction of the 
4d conformal plasma has equation of state defined by [21] 

»^4„-(^)». r.(^)'.„,3) 

with pq and a constants. This equation of state is valid in or out of equilibrium and is normalized 
such that the vacuum pressure vanishes. For the SS plasma in equilibrium, it follows from (|2.1ip 
and ()2.13p that the pressure and energy density satisfy the relations 

P = ^7^-po, P = P*7^ + Po, (2.14) 

where is a constant. 



2.3 Conserved charges 

The constituent fluid of the plasma object has local energy density p, pressure P, velocity n'^, local 
entropy density s, and local temperature T. These local quantities provide the information we need 
to compute the thermodynamic quantities (energy, angular momentum, entropy, temperature) of 
the plasma balls and plasma rings. 

To define these quantities recall that our fluid lives in a 3d Minkowski background with station- 
arity timelike Killing vector ^ = df and spacelike Killing vector x = d(f>- We can then foliate the 
spacetime into constant t hypersufaces and is their unit normal vector. Then, given any Killing 
vector if)^, one can define the associated conserved charges Q[^] = Jj^dV T^^C^ip'^, where dV is the 
induced volume measure on 'Et- The fluid velocity is given by ()2.10p . i.e., = ^ [6^^ + i^d'^^^ with 
7 = (1 -r^O^)-^/^ The energy and angular momentum of the plasma associated, respectively, 
with the Killing vectors ^ and x then 

E = [ dV[j^p + r^n^P)+a\df\5if)] , 
Jv 

J= [ dVr^n^^{p + P)-a [ dV h^.exldmf) . (2.15) 

Jv J'Et 

Note that for axisymmetric solutions (plasma balls and plasma rings) the boundary term in J 
proportional to a vanishes. It is however present for non-axisymmetric solutions where x ' ^ ¥^ ^ 
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(i.e., when the fluid boundary is not invariant under the action of x) i28j. The total entropy of the 
fluid is the conserved charge associated to the entropy density current su^, 



S = - dV{x-u)s= / dV-fs . (2.16) 
Jv Jv 

3 Equilibrium solutions: rigidly rotating plasma balls and rings 

There are three families of axisymmetric rigidly rotating equilibrium configurations in a 3d Minkowski 
background: plasma balls, plasma rings and plasma tubes. The later were already analyzed in a 
previous paper and do not interest us here [28] . The plasma balls and rings were discussed in detail 
in [21]. Because we will later study the stability of these solutions we review were their properties, 
following closely [21]. 

Consider plasma configurations in a d = 3 Minkowski background parametrized by coordinates 
(t, r, (/)). The axisymmetry requirement demands that the boundaries of the plasma depend only 
on r. Each boundary is thus defined by the condition (j specifies a particular boundary in the case 
where more than one is present) 

f{r)=r-Rj = 0, (3.1) 

d f 1 
and has unit normal = = 5^r- Its extrinsic curvature K = h^Vyu'^ Is K = 

Following [21], it is convenient to frame our discussion in terms of the dimensionless variables, 

^=^^ f=P^^ y = ftr = nf, (3.2) 
Po o- 

and also to use dimensionless thermodynamic quantities, 

E^esE. j^eM. s = -^^. f=T(°)"\ (3.3) 

We now consider the properties of plasma balls and rings. 
3.1 Plasma Balls 

Plasma balls are characterized by having a single axisymmetric outer surface at r = Rq and by 
P> = 0. Using equation of state (j2.13p . the Young-Laplace boundary condition (|2.7p reads 

piR,) = 4po + ^ . (3.4) 

Plasma balls in equilibrium must satisfy the equation of state (|2.14p and obey the boundary con- 
dition (j3.4p . This implies that 
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The range of v is [0, 1] and p{v) — po is always positive for the plasma ball. This is in agreement 
with the requirement that the local temperature defined in p.l3p must be positive. After using 
the equation of state (|2.13p and (|3.5p we get for the local plasma temperature, 

T = j[^P^^^y\ (3.6) 

For later use note that the constant introduced in (j2.14p is related to po by 

P* = ^Po9+(vo) ■ (3.7) 

This relation follows from replacing (12.14p in the Young-Laplace equation. 

Use of ()3.5p and (j2.13p in (|2.15p - (j3.3p yields for the dimensionless energy, angular momentum 
and entropy of the plasma ball, 

~^4jl-vt±^v^ J=M±^, S='SVY-^Jl + ^Y\ (3.8) 
while the dimensionless temperature and dimensionless angular velocity of the plasma balls are 

Note that the plasma ball temperature is the redshifted local temperature, T = T/7, in agreement 
with the discussion associated with (j2.1ip . The plasma angular velocity is naturally the same as 
the fluid one with no associated Lorentz factor. 

3.2 Plasma Rings 

These have an axisymmetric inner surface at r = R[ (where P< = 0), in addition to the outer 
surface at r = Ro (where P> =0). Using equation of state (|2.13p . the Young-Laplace equation 
yields the pair of boundary conditions, 

p{Ro) = 4po + ^ , p{Ri) = 4/50 - ^ . (3.10) 

Plasma rings in equilibrium must also satisfy the equation of state (I2.14p and obey the boundary 
conditions (|3.1Up . This means that rings must satisfy the pair of equations 

Piv)-P0/^ 2\2 (i , ^ Wi 2n2 _ , . 

[l-v) = [l + —\(l-v^] =g+{vo) 



3/90 V ^^o 



{l-vff^g.{vO. (3.11) 



V 



Note that p{v) — po, and thus the local temperature are non-negative as long as Vi > This system 
can be satisfied only when 

g+{vo) = g^Vi) . (3.12) 
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This condition constrains the three variables Vq, Vi and fl as, e.g., Vi = Vi{vo,Q). An inspection of 
g+{vo) = g-{vi) concludes that there is a minimum Vo, call it v*, above which ()3.12p is valid [2T] . 
So, plasma rings exist only for Vq > v*. In fact there are two families of black rings. One is called 
the fat plasma ring and exists for ^ < Vi < v* (where v* < v* is such that the derivative of g_(fi) 
vanishes), while the second, dubbed as thin plasma ring, exists for v* < Vi < 1. At Vi = v* the two 
families meet at a regular solution. 

Use of (j2.13p and (j3.12p yields for the local plasma temperature 



T 



\ a J \ a 



(3.13) 



Finally note that defined in (j2.14p is related to po by = 3pog+{vo) = 3pog-{vi). 

Using (j3.1ip and (j2.13p in (j2.15p - (|3.3p yields for the dimensionless energy, angular momentum 
and entropy 



E 



S 



+ 5n{vo + Vi) - n{v^ + vf 



J 



2{vt-vt) + 2n{v^ + vf) 



4 

^2 




f]3 



3/4' 



(3.14) 



while the dimensionless temperature and dimensionless angular velocity of the plasma rings are 



[g+ivo)]'/^ = [g^ivi)]'/\ 




(3.15) 



Note that to determine the plasma ring family of solutions, the constraint (j3.12p must be imposed. 
The different equilibrium configurations are best understood by looking at a phase diagram of 
solutions shown in later sections. 



3.3 The hydrodynamic regime 

Relativistic hydrodynamics provides a good effective description of the deconfined plasma phase of 
M = 4 Yang Mills theory compactified down to d = 3 on a Scherk-Schwarz circle only if certain 
conditions are satisfied [21]. First, hydrodynamics is by definition valid when the thermodynamic 
quantities of the fluid vary slowly over the mean free path i^jp of the fluid. In our case £mfp ~ ~ 
^ . A good estimate for the validity regime is obtained when the maximum fractional rate of change 
of the fluid local temperature, ^jj^^g^^ ~ dr IiitIuj^x (^^"^^^^ that T = T7) is much smaller than £mfp- 
This occurs for <C 1. This condition is satisfied by a wide range of plasma balls and rings 
as long as we are away from extremality. Second, the analysis done so far and onwards assumes 
that surface tension is constant, a = (7{Tc), when in fact it is a function of the fluid temperature 
at the surface. This assumption is valid when T/Tc ~ 1 at the boundary surfaces. This is the 
case for a long range of energies and angular momentum as long as we do not approach too much 
the extremal configurations. Finally, the boundary of the plasma is treated as a delta-like surface 
when in fact it has a thickness of order T^^. So the analysis is valid when the boundary radius is 
everywhere large when compared with T^^, {Rq, Ri, Ro — Ri} ^ This is the case if the plasma 
energy is large and again if we are away from the extremal configurations |21] . 
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4 Stability analysis of plasma balls and rings. 



In this section we want to consider a rigidly rotating plasma ball in d = 3 and address its stability 
when perturbed. The dynamics of the perturbations is dictated by the hydrodynamic equations, 
subject to appropriate boundary conditions. In this section we restore light velocity factors to ease 
comparison between relativistic and classical results. 

4.1 Perturbations of the equilibrium solutions: inviscid relativistic case 

Perturbations take the plasma away from thermal equilibrium and therefore viscosity and diffusion 
effects start to contribute. The energy-momentum tensor of the fluid includes not only the perfect 
fluid and the boundary surface tension terms (j2.ip . but also a dissipative contribution (j2.2p . For 
now, we neglect the dissipation contribution to the fluid stress tensor. 

Consider then a generic equilibrium solution described by velocity u'^q^ = 7 [6^^ + ^6'^'^) (with 

pressure and density functions as given by ()2.14p . which we label with 
the subscript (0), standing for unperturbed quantities. Now, suppose the system acted upon by a 
perturbation with the generic form 

P = P(o) + 6P , 6P{t, r, (j)) = ^^^■p(r)e^-^-ima)t+im^ ^ ^4 -^^ 

3 3 

P = P(0) + , P(Q) = P{0) + ^Po , Sp{t, r, (/)) = -2 SP{t, r, 4>) , 

where we used the equation of state (j2.13p valid also out of equilibrium and we denote the pertur- 
bation of a quantity Q as 6Q. Positive real part oiu signals an instability. After eliminating the 0*^ 
order terms using the unperturbed hydrodynamic equations, the continuity and the Navier-Stokes 
equations yield, up to first order in the perturbation, 

c^nJj^V^^p + c^6u''Vf,p^o) + (P(o)C^ + P{o))V^,Su^' + {c^6p + 5P)V^ufo) = , 



+ (g^" + <o)<o)) Vm'^P + (Su-'u'l^^ + u'(,^6u^') V^P(o) = . (4.2) 



which reads 



3 Tuj d 

^'^^ V + -f'^ — {rUr) +r{uj - imVL)Ut + imU^ , (4.3) 



' dr 



^ [iujVLr'^ - mc'^-f-^)V - 2n^r-f-^Ur - c^co Ut , (4.4) 



4p» „2 
3 



= T^V +ujUr-2^''m^, (4.5) 
3 

= -^^[iooVlr'^ -mc^-f-'^)V -2Q.r-f-^Ur -ivrU^, (4.6) 
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where V' = Multiplying (USD by Q and subtracting (031) it follows that Ut = '-§-U^, which 
satisfies the requirement that u^n^ = — up to order e. These equations can be used to get a 
second order ODE for V, and another equation defining Ur in terms of V and its derivative. The 
later is 

-2in \c^m - r'^nUuJ + ml7)l -f'^'P(r) - c^ujrV'(r) 
UJr) = ^ — — , (4.7 

while the second order ODE for V is 

= - cS"^ (c2 u;2^-2 + 4^2^72^ p// _ ^^6^-4 (^2(^2 ^ j^2 ^2) ^ 4^2j^2) 

+ c^m2(Lj2 + 41^2) _ ^;2 ^8j^6(-^ _ -^^^^2 _^ g4^4j^2 (|g^^3 ^ _ _^ 2u;2f]2(3m2 - 14) + 24iu; 



+ c^r^ {Suj'^ - 2iu;^ mCl- Auj^{m^ - 5)1^^ - 12iiu; mQ^ + 4(8 - 3m^)n'^) 



V 



The perturbed continuity and Navier-Stokes equations must be supplemented by appropriate 
boundary conditions. For that, let us write a general boundary disturbance as 

r = R{t, (t>) , with R{t, (t>) = Rj{l + ex e(^-^™f^)*+^™<^) , e < 1 , (4.9) 

where Rj is the unperturbed radius of the boundary. For the plasma ball one has Rj = Rq, while 
for the plasma ring we have both the inner and outer boundaries: Rj = R^ and Rj = Rq- For 
plasma rings, this reduces the possible set of perturbations to the subsector with similar temporal 
and angular deformations in both boundaries. 

The first boundary condition is a kinematic condition requiring that the normal component of 
the fluid velocity on the boundary satisfles the perturbed version of (|2.8|) . u^^-^ 5n^ + 5u^rS^^ = 0, 

where 5n^ = — nf^^ and the unperturbed normal is n'^^ = ~ "^a*' 'This ensures that 

the fluid is confined inside the boundary and is also a consistency relation between the boundary 
and the velocity perturbation. To leading order this boundary condition reads, 

BCI: Ur\^^.^uj^t,iRo^Xo,i. (4.10) 

The second boundary condition is a balance on the normal stress at the boundary. This means that 
the pressure perturbation must also satisfy the perturbed version of the Young-Laplace equation 
(|2.7|) . namely: (5P< — 5Py = a5K. Since we have vacuum in the exterior of the plasma configuration 
this reads 

BCII: iP^<^\uu...+5P\u ^-P<^\r> =±a(K\^,,^,-K\ 



$ \R{t,<t>)^ \Ros) \Ro,i V '^(*'<A) ^Ro,i, 

where the choices {<,+,o} apply to the outer boundary and {>,— ,i} to the inner boundary, if 
present. The subscript r = R{t, (p) means that we evaluate the expression at the perturbed bound- 
ary r = R{t, (j)) defined in (j4.9p and the subscript r = Ro^i means evaluation at the unperturbed 
boundary r = Ro,\- is computed using (j2.14p . The extrinsic curvature K = h^'S/,^n'^ is 

obtained using the unit normal of (j4.9p . 

n, = 15/1-1 + 5; - R'^5^) , \df\ = (^1 - + }_j^>2^ ' (411) 
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The Young-Laplace equation ()4.10p then yields to leading order in e 



-nl ^ -4 



(4.12) 



where it applies, respectively, to the outer and inner boundary and 7oi = tIo • Here, we have 
defined the rotational Bond parameter S which plays an important role in this problem. It measures 
the competition between centrifugal and surface tension effects and is defined by, 

S = ^Rl, ^ . (4.13) 

4.2 Plasma balls: instability and critical rotation in the inviscid case 

We now particularize the above framework for plasma balls. We can combine boundary conditions 
(j4.12p and (j4.10p in a single condition, 

UriRo) = ^ V{Ro) ■ (4.14) 

a I — 1 — H — t{u! — imO) 



? 

To summarize, (j4.7p and (j4.14p give us a condition on V. Together with regularity conditions at 
the origin, Eq. (j4.8p is then an eigenvalue problem for uj. 

Although we will present our numerical results in full generality, it is insightful to compare 
them with the small rotation regime where an analytical treatment is possible. In this small 
velocity regime, Q.R0 <C c, (|4.8p . (|4.7p and ()4.14p reduce respectively to 

J|:(-n-^^ = o, (4.15) 
^.(r) = -^T'^^^^^"^^^^^\ (4.16) 

UriRo) = , V{Ro) . (4.17) 

a (m^ — 1 — Zj) 



If we define 



T^^^^RW. (4.18) 



we can express the analytical solution to the system ()4.15p - ()4.17p as 

^/tZ = iVr, ± y^{m - 1) (S - m{m + 1)) . (4.19) 



Thus, for S > m(m + 1) and m > 2, the system is unstable. This is one of our main results: in 
the simplest m = 2 case, plasma balls become unstable against 2-lobed or peanut-like deformations 
when the rotation reaches a critical value Ocrit- For higher rotation they become also unstable 
against m-lobed deformations, with m> 2. Notice that in the non-relativistic regime the density is 
approximately a constant and equal to ^p^,/2> (see footnote [5]) . Thus, our result ()4.19p is precisely 
the well-known result for inviscid, incompressible fluids of density p = 4p*/3 |31j . 
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Figure 3: Left Panel: Details of instability for the m = 2 mode. We are fixing E = I X 10^, 
a = 10"''. Classically, i.e., for small rotation rates, the quantity is a constant, and both the 
real and imaginary part of the characteristic frequency uj lie on the straight line, i.e., directly 
proportional to The stability region lies in a range of very small QR^, not visible in the figure. 
Right Panel: Details of the instability for the m = 2 mode, this time with parameters chosen such 
that the instability sets in at large 0,Ro. In this case, the threshold is around VlRo/c ~ 0.402, the 
parameters are = 24. Classically, the threshold would be at at VLRq/c = 0.5. 



Our numerical results for the relativistic system (14.81) . (14. 7p and (|4.14p . are depicted in Fig. [3] 
for the m = 2 case. For small rotation rates, they are in perfect agreement with the non-relativistic 
limit (j4.19p . We find a (in any case small) deviation from the classical prediction only when VtRo/c 
approaches unity. In the right panel, we show a case where the threshold rotation frequency at 
which an instability sets in is rather large. For these values we get a threshold of approximately 
VLRo/c ~ 0.402, still in good agreement with the classical result VLRq/c = 0.5, as given by (j4.19p 
for these values = 24). We were not able to find unstable modes for m = 1, in agreement 

with the classical result (|4.19p for small rotations. 

The critical rotation frequency i?of^crit) for which the configuration is marginally stable is shown 
in Figure m For a given p^Roja and rotations larger than Ro^ccit^ the system is unstable to two 
lobed perturbations [m = 2). We find that for moderately large p^Roja > 50, the classical formula 
(I419D holds. 

At this point we should check if and when the lengthscale of the instability falls within the 
hydrodynamic limit (discussed in the end of Section ED . The thermodynamic quantities of the fluid 
must vary slowly over the mean free path of our system, £jnfp ~ ~ The lenghscale of the 
instability is the ratio of its wavelength to the tube's radius and is of order (lj/?^)^^. Therefore the 
hydrodynamic description is valid for instabilities that satisfy ujRq 3> — ^r^- ^ simple inspection 
of Fig. [3] allows one to conclude that the conditions for a hydrodynamic description are satisfied in 
the non-relativistic limit where p*'^J^° > 50. It is certainly satisfied for very large values of IM-Es.^ 
In addition, the initially unperturbed plasma ball must of course satisfy the conditions discussed 
already in the end of Section [3l 
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Figure 4: Critical rotation frequency ^IRg at which instabihty sets in, as a function of the dimen- 
sionless quantity ^3^- For large ^^3^° , the classical prediction ()4.19p applies. 

4.3 Viscosity: instability and critical rotation in the non-relativistic case 

Although we start with a plasma ball in thermodynamic equilibrium, viscosity contributions have 
to be taken into account in the perturbed configuration. In fact even the most "perfect" known 
fluids have a certain non- vanishing viscosity [32], and this is certainly the case for the SS plasma. 
It turns out that, for the problem at hand, even a vanishingly small viscosity has a dramatic effect 
on the critical rotation at which instability sets in. In fact, viscosity introduces a singular-limit, 
where the "limit of theory is not the theory of the limit" (see [3H [33] for the classical analysis). 
Thus our relativistic analysis of the previous subsection must address also viscosity effects. 

Fortunately, from the full relativistic analysis of the previous section, we know that we can 
to a good approximation address the present problem in the small rotation regime. To define 
precisely this regime, recall that in section 13.31 we concluded that hydrodynamics provides a good 
effective description of the deconfined plasma phase immersed in the vacuum confined phase when 
the plasma satisfies the condition 

«1. (4.20) 

Moreover, in the previous subsection we found that the full (numerical) results for the marginally 
stable mode (where the m-lobed instability set in) agree very well with the (analytical) non- 
relativistic results obtained for 

S IRn , N 

-<^1. (4.21) 

c 

We are therefore justified to use the small velocity regime: in this case not only the non-relativistic 
results reveal very good agreement with the full relativist analysis but also this is the relevant 
regime where the hydrodynamic analysis provides valuable information for the dual gravitational 
system. 

The non-relativistic limit of the hydrodynamic equations presented in subsection 12.11 was studied 
in great detail in [34] . There it is found that the continuity and Navier-Stokes equations reduce. 
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respectively, to 



V • V = 0, 

dtv + ( V V) V = - — VP + lyV^v , 1^=^, (4.22) 
4/0* 4p^, 

where one uses u — > (l,v), and Vj (i = r, (/>) represents the covariant derivative with respect to 
the purely spatial metric ijij. So in the non-relativistic limit, the hydrodynamic system reduces to 
the continuity and Navier-Stokes equations for an incompressible fluid with constant density p\^_i 
(so the continuity equation simply states that the velocity field is a solenoid vector: its spatial 
divergence vanishes), and with kinematical viscosity u^. 

At this point we can now perturb (|4.22[) and simply follow the classical analysis of the m- 
lobed instability of a fluid ball done e.g., in [SJj. Using the same non-relativistic version of the 
perturbations used in (|4.2p . we get. 



-^(r Ur) + imllA, = , 
ar 

A!^p + 217f/. + .f/, = .(f/; + if/;-(!^^[/, + ^f/.) . (4.23) 

These equations must now be supplemented by the appropriate boundary conditions. The con- 
straint ()2.8p is still valid as well as the associated boundary condition (14.10(1 . The Young-Laplace 
equation gets now a contribution from the viscous term, and so the normal stress-balance at the 
boundary, ()4.12p . is modified to 



[m 



2 1^ 

^ 3 a 



+ 2vU'^. (4.24) 



Finally, we must also require that the tangential stresses vanish at the boundary (this amounts to 
require that the fiuid is shearless at the boundary), which yields the extra boundary condition 

RlU'^ - RoU^ - imUr = . (4.25) 

We can use the first relation in (|4.23p to express Uff, in terms of Ur and its derivative and then use 
the remaining equations to solve for Ur and V. One gets a fourth-order differential equation. In 
any case, the procedure is trivial and one ends up with the following eigenvalue equation [3T| 



l3^ + 2l3'^(^{rr? -l)-m{m-lf^-iR(^ +mRe^ (^ ""^^ ^ - 1^ =0, (4.26) 



It follows from (|3.5|) and (|3.7p that in the non-relativistic limit one has p ~ |p, ~ 4po- Moreover, in this limit 
one finds that a possible contribution coming from the bulk viscosity in the Navier-Stokes vanishes because it is 
proportional to V ■ v, and the absence of the particle number conservation and the use of the Landau frame [ID] allow 
to avoid the use of thermal conductivity [34] , 
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where we defined 



Re^^, P^^^Re, $ ^ ^ 2m - 2/3:^!Ii^iM±#±lM . (4.27) 

Here Im{P) is a modified Bessel function and Re is the Reynolds number. Solving for this eigenvalue 
equation, we get the behavior depicted in Fig. [5l What the figure shows, and can be proven be 




Figure 5: Instability as a function of the inverse of the Bond number, for different values of the 
Reynolds number Re = 10, 100 and for m = 2. We also show the inviscid values. Notice that the 
critical Bond number is not m{m + 1) when viscosity is present, even in the limit where it vanishes. 
In the general viscous case, the instability is limited by the critical point Sc = — 1. 

proven analytically [31 j, is that the limit of small viscosity. Re oo, does not yield the zero 
viscosity result. This is quite an astonishing result: the critical Bond number Sc when viscosity is 
introduced is at Sc = m'^ — 1 that is always smaller than the critical value for the instability found 
in the last subsection, namely, m{m + 1). Therefore, the configuration is unstable at lower rotation 
frequencies, when viscosity is accounted for. This result holds for both large and small Reynolds 
number. 

In Fig. [6] we represent the plasma rings, and stable and unstable plasma balls in two distinct 
phase diagrams at fixed energjo, and it summarizes one of our main results. 



5 Bifurcation to two-lobed configurations: rotating plasma peanuts 

In the previous section we found that plasma balls are marginally stable at Sc = — 1, when 
perturbed by an azimuthal mode m. Such a marginal mode usually signals a branching off to 
another family of solutions, and here it is no exception. This new family is a non-axisymmetric 

^In these diagrams we choose to fix the dimensionless energy at the value i5 = 40 to make a connection with the 
value chosen in [21], where plasma balls and rings where first discussed. This value of the energy corresponds to 
— - ~ 14. In this case the non-relativistic analysis for the marginal point mismatches the accurate full relativistic 
result by a factor of approximately 20%. 
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Figure 6: (a) Phase diagram with the entropy of plasma balls and plasma rings E as a function 
of angular momentum J, at fixed energy E = 40 in d = 3 [21] • For this value of the energy, one 
has V* ~ 0.7248999 and v* ~ 0.5009615. Plasma rings exist only for Vq > v*. The fat plasma 
rings exist for Q{vo = 1) < i"! < v* (with i}{vo = 1) = 0.333333), while the thin plasma rings 
exist for v* < Vi < 1. At Vi = v* the two families meet in a regular solution. Hydrodynamics does 
not provide a good dual description for configurations near extremality (T = 0) where the entropy 
also vanishes. We find that plasma balls become unstable above the critical rotation Jc — 21.8127, 
and the dashed line represents the unstable plasma balls, (b) Similar to Fig a), but this time we 
represent the phase diagram with the angular velocity of plasma balls and plasma rings as a 
function of angular momentum J, at fixed energy, E = 40. 

m— lobed configuration: a rotating plasma peanut. In this section we want to verify the existence of 
this branch of solutions and study some of its main properties. A detailed numerical study of the 
full branch of this new plasma phase is outside of the main scope of this work. Here, we will address 
the most important region, namely we investigate how this family looks close to the bifurcation 
point S = Sc in the phase diagram of stationary solutions. Consistent with the analysis of previous 
sections we restrict our analysis to the hydrodynamic and small rotation regime, Eqs. (|4.20p and 
(I4.2ip . where the hydrodynamic analysis provides valuable information for the dual gravitational 
system. We can then follow closely the analysis done by Benner, Basaran and Scriven [33]. In the 
following we obtain the new branch of rotating plasma peanuts that emerges from the plasma ball 
bifurcation point. 

The energy, angular momentum and entropy of the plasma lump are given by (j2.15p and (j2.16p . 
In the regime (j4.20p and (j4.2ip (see also footnote [5] and note that pc^ + P ~ 4poc^) and in the close 
vicinity of the plasma ball bifurcation point the associated dimensionless charges read, 

^ = ^ / T« « 1 r d<t>R\e,ct>), 

T^O Jv ^ Jo 

^ = TTTT / 75 « - / d(t>R (e,(/)), 

J = [ r'^T"*' « - r d<pnie)R\e,<P). (5.1) 

TTo-^ Jv vr Jq 

In these relations we are looking for the unknown family of lobed solutions by expanding around the 
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known plasma ball at the marginal stability point Sc = — 1. We take e to be the parameter that 
measures the deviation from the rotating axisymmetric plasma ball with dimensionless radius Rq, 
and i?(e, (f)) describes the non-axisymmetric boundary of the unknown plasma peanut (note that in 
the vicinity of the bifurcation point the axisymmetric deviation is really small and the perturbative 
approach is appropriate). We also assume that there is mirror symmetry around the </> = 0,7r axis. 

In a perturbative analysis, known as power series method [33|, we now take the boundary radius 
and Bond parameter to be described by the expansion around the marginally stable plasma ball, 

R{eA) = i?o(l + e/'H<A) + f/')(0)) +0(e'), 

S(e) = Se + eS« + +0(e3). (5.2) 

We also introduce the dimensionless pressure jump IT at the axis of rotation (see ()2.14p ). and its 
expansion, 

n ^ ^ Q - po) , n(e) = n, + en^) + 1 n(2) + o{e^) . (5.3) 

The perturbed quantities characterize the m-lobed branch. We want to find a family of configu- 
rations that have fixed energy. Perturbations have to solve a total of four equations, namely: the 
Young-Laplace equation, the energy constraint (we want to represent the branch of solutions in a 
phase diagram at fixed energy) , the orthogonality equations and a equation defining the amplitude 
parameter e. 

The perturbed n^^ order Young-Laplace equation, 5P< = adK, yields 

/^) + (i + s,)/W + n(") = 5K, (5.4) 

where = d^f and the relevant source terms Syi are 

= 0, 4i = 0, 41 = (2-Se)(/«)'-2S«/W + (/«)V4/«/«, 
SP, = -3S(2)/(i)+3(2-Se)/«/(2)+3(/«/f -6(/«)VW 

+6 ( f^'^ffl + /(^)/«) + 9/«(/r )^ - 6(/f)3 - 9/(^)(/«)^ . (5.5) 



The condition that fixes the energy follows from perturbation of the first relation in (j5.ip . 5E = 0, 
which yields 

rd^f(-) = SP, (5.6) 
Jo 

with source terms 

4°^=0, 5«=0, = - [f^'))\ sP=3l\<PfWf(^K (5.7) 

For n > 2 the problem is inhomogeneous and the solution must satisfy an orthogonality conditioE0, 

r #/(^)5S + 2n(i)4") = . (5.8) 
Jo 

^To get (|5.8p start with the Young-Laplace equation H5.4|) . Multiply it by Z*-^' (</>); integrate over <j!>; do a double 
integration by parts (use the symmetry condition /'"'(O) = /''"■' (tt) and the n = 1 Young-Laplace to simplify some 
of the terms); and finally make use of the energy conservation (|5.6p . 
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Finally, the amplitude parameter e is defined as the integral-weighted difference between plasma 
shapes, 

e^£d<PfW(^R-R,) ^ l\cPfWf(^) = 6i. (5.9) 
The solution of (|5.4p - (l5.9p . up to second order in the perturbation, is 



/«(0) = 


Y I cos(mc 
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^ m? — I 
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^ — 1 ) cos^(m0) H — ( — ^ — 3 ) sin^(m(/)) , 
71^ J IT J 

m'^ - 1)(1 -m^) 

(1 - m2) . (5.10) 

^ z vr 

These perturbations keep the energy of the stationary solutions fixed and equal to E' = 4i?^. 
Moreover, use of these relations in the definition (I4.13|) of the Bond parameter and in (15. ip yields 
the expansion of the other thermodynamic quantities around the bifurcation point. 



S = ARl + O (e^) . (5.11) 

This expansion allows to represent the new branch of non-axisymmetric plasma peanuts in a phase 
diagram of stationary solutions at fixed energy. Since rotation is the mechanism responsible for the 
instability that signals the bifurcation to the new phase, it is appropriate to represent the 2-lobed 
plasma lump in a phase diagram that represents the angular momentum J against its conjugated 
chemical potential Vt. To check the accuracy of our approximations we first observe that at zero 
order, we have J^^^ = 2R'^0, = j^E'^Q. For the slope of the plasma ball at the bifurcation point we 
thus get gg- = for ^ = 40. This value is in reasonable agreement with the exact relativistic 
value for the $7 — J slope at the bifurcation, obtained from (j3.8p (for ii^ = 40 we get numerically 
. . 1.2 . explained in footnote [6] the disagreement is due to the fact that the system with 



dj(0) 200 ' 

£^ = 40 is not exactly within the classical regime). We can now use the next-to-leading order 
non-vanishing contribution to get the desired slope for the m-lobed branch, yielding 

2A{m'^ - 1) 1 




3 + 32m2 - 3m4 £;2 ' 



(5.12) 



For m = 2 this is the slope at which the 2-lobed or plasma peanut branch emerges from the 
plasma ball bifurcation point. In the phase diagram $7 — J of stationary solutions at fixed energy, 
represented in Fig. [71 the plasma peanut branch bends down and to the right relatively to the 
plasma ball family. In this diagram we identify the critical unstable point where bifurcation occurs 
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and we represent by a dashed line the plasma balls rotating faster than the critical velocity, which 
are therefore unstable. 

This figure summarizes the main result of this section: we have confirmed the existence of a 
new phase of non-axisymmetric stationary solutions and we were able to use perturbative methods 
around the bifurcation point to find what is the direction that the new branch of solutions takes 
relatively to the known phases. A full description of the 2-lobed branch well away from the bifurca- 
tion point (where it acquires a well defined peanut shape) would required a full numerical analysis 
and we leave it for future work. Note also that keeping going up along the plasma ball branch this 
time already in the unstable ball region to 2-lobed perturbations we would find a succession of new 
bifurcation points to new phases of solutions representing m-lobed plasma lumps with m > 2. The 
slope of these branches is given by ()5.12p . 




Figure 7: Phase diagram 0,(J) of stationary plasma solutions with fixed energy, E = 40. At 
(Jc^c) — (21.81,0.157) the plasma ball becomes unstable and bifurcates to a branch of rotating 
plasma peanuts. We show this branch in the vicinity of the bifurcation point but it continues for 
larger values of J. Note that the bifurcation point for plasma balls is at lower angular velocity and 
angular momentum than the merging point between the fat and thin plasma rings, a feature we 
verified to be valid independently of the choice for the energy of the system. 

Alternatively, we could represent the new phase in the S—J phase diagram. We get {dS/dJY'^^ = 
up to second order terms (this is not a surprise since at the order we work, the entropy is equal 
to the energy). Now, in the classical regime, the plasma family also has zero slope for this quantity 
so in principle one should go to higher order in e if we wish to compute accurately the slope of the 
m— lobed family. In our approximation the bifurcation point tends to collapse to the static plasma 
ball point in the S — J diagram where the slope vanishes. Therefore, at the order we work, the 
representation in this diagram is not a good choice. 

The current line of research is not exhausted in the present analysis. In future work, it should 
be certainly possible to find numerically the full branch of non-axisymmetric plasma lumps in the 
phase diagram, and to discuss their stability. 
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6 Unstable plasma balls and their dual black holes 



In the previous two sections we found that, starting with a static plasma ball, as we increase its 
rotation, the axisymmetric rotating plasma ball becomes unstable first to a 2-lobed perturbation 
and then to m > 2 lobed perturbations. Moreover, the marginal unstable modes were found to be 
bifurcation points in the phase diagram of solutions to a new branch of m-lobed plasma lumps. It 
is important to emphasize that these new plasma phases inhabit a region in parameter space where 
hydrodynamics provides a reliable holographic description of the dual gravitational system. This 
dual system was discussed in the Introduction: in the long wavelength regime, 3d fluid dynamics is 
an effective theory describing the Scherk-Schwarz (SS) compactification of a 4d CFT. The latter is 
dual to the SS compactification of AdSs gravity. It then follows that 3d plasma lumps correspond, 
in the dual gravity description, to SS AdSs black objects. The axisymetric plasma balls and plasma 
rings correspond, respectively, to rotating black holes and black rings in SS AdS5 [21]. 

Since plasma balls are unstable above a critical rotation rate (cf. Section^), the holographic dual 
SS AdSs black holes must also be unstable against m-lobed perturbations. These axisymmetric 
black holes are expected to bifurcate to a new branch of non-axisymmetric black hole solutions 
in accordance with the plasma results. In the simplest m = 2 case, these are non-axisymmetric 
solutions that describe a peanut-like rigidly rotating black hole. Such a deformed black hole must 
necessarily emit all kinds of radiation. These waves can escape to infinity through the directions 
parallel to the holographic boundary (where the dual fluid lives). It is then natural to expect 
that this non-axisymmetric black hole will decay, probably emitting a good amount of the angular 
momentum and multipoles of the system, into an axisymmetric slowly rotating black holeo In 
practice, this instability provides then a mechanism that constraints the rotation of the black hole: 
effectively it introduces an upper bound for the rotation of the SS black hole. Some of these 
non-axisymmetric black holes are expected to decay very slowly and to be long-lived. Indeed it is 
important to emphasize that hydrodynamics can only provide a good description of a gravitational 
system in a regime where the gravitational interaction and radiation is suppressed. For example, 
two plasma balls do not interact and their collision is not accompanied by radiation emission, so 
they can only describe approximately two black holes when they are widely separated. Therefore, in 
the regime where our fluid description provides a good approximation, we should expect radiation 
emission to be suppressed in the dual gravitational system. Thus the non-axisymmetric black holes 
should indeed be long-lived, leaking very slowly radiation and angular momentum. Finding the 
explicit geometry describing such a solution will probably require a full-blown numerical solution 
of the field equations. This is thus a valuable example of the power of the hydrodynamic/gravity 
duality. 

We expect our main conclusions to extend to higher dimensional theories as well. For any 
d > 3, (i-dimensional fluid dynamics is an effective theory describing the SS compactification of 
a {d + l)-dimensional CFT, which is dual to a SS compactification of AdSrf-|-2- Apart from an 
expected simplification of the technical analysis in the d = 3 case, there is clearly no step in our 
analysis that is valid only for d = 3. Thus, plasma balls and dual black holes rotating above a 

*We are undoubtedly grateful to Roberto Emparan for fruitful discussions concerning the dual interpretation of 
the plasma results 
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critical rate should also be unstable in higher dimensions and bifurcate to m-lobed plasmas and 
black holes. But, as we discuss in the sequel, we also expect some significant differences between 
the d = 3 and d > 3 cases. 

Eventually new features are expected as one climbs the dimensionality ladder. For instance, 
new plasma lump configurations appear for d > 3. Indeed, in d > 3 one can have, besides the con- 
figurations discussed here, also pinched plasma balls, non-uniform tubes and pinched non-uniform 
tubes [211 [22l [28] . In the dual gravitational system, pinched balls [2^ correspond to pinched black 
holes |35j . and the non-uniform tubes describe non-uniform black strings with the pinched ones 
being black strings on the verge of expelling a black ring [28]. It would be quite interesting to 
understand the stability properties and phase diagrams of plasma balls in higher dimensions. 

But, the results of the present analysis already allow one to infer important properties about the 
stability and existence of new solutions in higher dimensions. Consider hydrodynamics in a d = 4 
Minkowski background (the boundary of SS AdSe black objects). Take a static plasma tube with 
topology X M or x (it is translationally invariant along the extra direction which can be 
or not compact). Such a plasma tube is unstable against the Rayleigh-Plateau instability when its 
length is larger than its transverse perimeter [HO [28]. In the dual system, the corresponding black 
string is Gregory-Laflamme unstable [3] , the holographic dual of the fluid instability. The marginal 
unstable mode is a bifurcation point to a new phase of solutions describing non-uniform plasma 
tubes and black strings. For details of this instability, see Refs. [11 [Sj [HI [28l [36] . Now, a Ad static 
plasma tube is simply a 3d plasma ball trivially extended along the extra direction. Our results 
then show that a rotating plasma tube, with the rotation axis along the tube direction, should also 
be unstable, above a critical rotation rate, to m— lobed perturbations. Moreover, the marginally 
stable points are bifurcation points, this time to non-axisymmetric plasma tubes translationally 
invariant along the tube direction and whose transverse cross section has a m-lobed shape. Again, 
by the hydrodynamic/gravity duality, we expect that rotating black strings will become unstable, 
not only against the Gregory-Laflamme instability, but also against m-lobed azimuthal perturba- 
tions. And a new branch of non-axisymmetric m-lobed black strings is expected to branch-off at the 
unstable threshold point in the phase diagram of solutions. Other interesting solutions include ax- 
isymmetric non-uniform black strings (these bifurcate from the Gregory-Laflamme unstable point) 
and pinched non- uniform black strings which are also solutions of the gravity theory [28]. Our 
study then predicts the existence of non-axisymmetric non-uniform black strings and pinched ones. 
We emphasize that the existence of pinched black strings [28j and non-axisymmetric m-lobed black 
strings would hardly be anticipated without resorting to the hydrodynamic/gravity duality. Our 
discussion focused on 4d plasma tubes dual to 6d black strings but again it should extend to higher 
dimensions. 

These results rely on the hydrodynamic description of a particular gauge/gravity duality, namely 
of the Scherk-Schwarz system. Experience with hydrodynamics indicates that similar results should 
be found for other gauge/gravity theories with an effective hydrodynamic description and a confine- 
ment/deconfinement phase transition. The deconfined black hole phase of such a general system is 
still expected to be described at leading order by a perfect fluid holographic stress tensor, and the 
interface between the deconfined and confined phases is again expected to be dictated by a domain 
wall with a surface tension. Such a generic dual system will of course have a different equation 
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of state, encoding the information on the kind of fluid describing the gauge theory. However, our 
analysis is not very sensitive to the particular equation of state of the plasma. Hence it could be 
that our main results on the stability and bifurcation properties of plasma balls will be common to 
generic dualities. 

It would also be interesting to investigate the dual of the superradiant instability on the plasma 
balls. In general, rotating black holes develop an ergoregion, a region in spacetime with negative 
energy states. For spacetimes with ergoregions, one can have superradiant scattering, whereby a 
wave (with frequency to < mil, where m is an azimuthal number and O is the horizon velocity) 
can be amplified, extracting rotational energy from the hole. In AdS, the superradiantly amplified 
waves bounce back at infinity and lead to a superradiant instability [33]. Thus four- [421 [Hj and 
higher dimensional [45 1 [36| rotating Kerr- AdS black holes can be unstable. The endpoint of such 
an instability is presumably an element of a new branch of stationary black holes, rotating in such 
a way as to avoid the superradiant window [451 132] ■ Some of the angular momentum of the black 
hole is transferred to caged radiation in between the horizon and the AdS wall and co-rotating 
with the black hole. It was also argued in [45J that a new branch of non-axisymmetric black hole 
solutions could eventually bifurcate from the original Kerr- AdS black hole at the threshold of the 
superradiant instability. This superradiant phenomena is a similar, but not identical, mechanism 
to the one we explored in this paper (the mechanism dealt with here is not dependent on super- 
radiant amplification). They may correspond to different bifurcation branches in a phase diagram 
of possible solutions. To explore the eventual dual of the superradiant instability on a plasma ball 
we should take into account that experience indicates that an ergoregion instability develops when 
the rotation speed at the boundary surface, il.Ro, exceeds the sound speed of the plasma. This 
then corresponds to the formation of an acoustic "ergoregion" in the system, and therefore by a 
general theorem by Friedman |i47j, they should be unstable. We should however keep in mind a 
possible serious caveat: the dual hydrodynamic description is typically valid for large AdS black 
holes (i.e., for those whose horizon is large compared with the cosmological scale) |2H I30|. while 
the superradiant instability is present only in small AdS black holes [42^ I44j. So it might well be 
that the fluid description is not able to capture the dual of the superradiant instability. 
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Appendices 



A On the stability of plasma rings 

In this Appendix we discuss the stability of plasma rings, and we find that plasma rings are stable 
against the particular m-lobed deformations that we consider. Plasma rings have an inner boundary 
at r = i?i in addition to the outer boundary at r = Rq- The boundary conditions ()4.12p and (|4.10p 
yield the following system for the outer and inner surfaces, 

"^-"-0,1 1 0,1 



Ur{Ro,i) 



a 



± (^(u; - ^mnYRl, + - l) - Ifi ^ 76 . 



V{Ro,{) . (A.l) 



Relations (j4.7p and (lA.ip provide two conditions on V and constitute again an eigenvalue problem 
for uj. The problem depends on two dimensionless quantities: ^*^° and VlRo. For Q.R0 <C c, (j4.8p 
and ()4.7p reduce respectively to ()4.15p and (I4.16p . while (|A.ip simplifies to 

For black rings some of these quantities are not independent, e.g., we have i?i = R\{Ro) and 
Q, = il,{Ro, p/a). In section [3] we already discussed these constraints so it is appropriate to rewrite 
()4.15p , (j4.16p and (|A.2p in terms of the dimensionless variables ()3.2p . Defining also the dimensionless 
frequency uj = -^uj the system we have to solve is 

V{v) = Av'^ + Bv-"^ , (A.3) 
2imnV{vo)+uVoV'{vo) vl ^^^^^^ 



4t;5+(?;o) + 41^2^ [± (m2 - 1) - i5+(^^o)i^3 



4u5_(ui) (52 + 402 j n ± (m2 - 1) - lg_{vi)vl 



V{vi) . (A.5) 



We have found no unstable mode. One might argue that black rings should be unstable at 
least at the point where they cross that plasma ball diagram of Fig. [7l This does not necessarily 
have to be in contradiction with our results first because the boundary conditions for black rings 
do not seem to reduce to the boundary conditions for plasma balls, even in this limit. Second, 
the analysis here does not prove that these plasma rings are stable: they are stable against the 
particular m-lobed deformations that we consider (see discussion just after ()4.9p ). A more general 
analysis is needed, encompassing generic perturbations. 
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